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ABSTRACT 


The  current  etete-of-the-art  is  reviewed  with  respect  to  the  calculation  of 
convective  stagnation  point  heating  at  supersateUite  reentry  speeds.  Recently 
calculated  transport  properties  have  been  compare'^  with  experiment.  It  is 
noted  that  theoretical  and  experimental  eetimates  of  the  total  thermal  conductivity 
are  in  much  closer  agreement  than  reported  by  earlier  investigators. 

Similarity  solutions  employing  these  recently  computed  transport  properties 
are  presented  for  the  convective  heat  transfer  rate  in  an  ionized,  dissociated 
gas  for  equilibrium  air  and  equilibrium  nitrogen  at  reentry  speeds  up  to 
70, 000  feet/sec.  Solutions  are  also  obtained  for  the  case  when  large  rates  of 
injection  are  introduced  at  the  stagnation  point.  Tablet  of  boundary  layer 
cluracteristics  including  profiles  of  temperature,  velo:i'.ty,  and  enthalpy  are 
presented  for  the  axisymmetric  and  two>dimensional  stagnation  point. 

This  task  was  initiated  under  the  REST  Program  under  the  cognizance  of  Air 
Force  Ballistic  Systems  Division;  the  major  portion  of  the  work,  however, 
has  been  company  supported. 
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I.  INTRODUCTION 


One  of  the  immediate  problems  in  the  design  of  v'^Hicles  entering  the  Earth's 
atmosphere  at  supersatellite  velocities  is  the  determination  of  heat  transfer 
from  a  stream  containing  dissociated  and  partially  ionized  gas.  In  order  to 
attack  this  problem,  one  must  know  with  a  certain  degree  of  accuracy  the  trans¬ 
port  and  thermodynamic  properties,  and  the  rate  processes  involved.  A  complete 
solution  of  the  overall  problem,  where  one  takes  into  account  body  geometry  and 
general  ablating  surfaces,  is  at  present  nonexistent.  Specialized  solutions  are 
available,  however.  Examples  are:  simplified  chemistry  (i. e. ,  equilibrium), 
simplified  geometry  (i.e. ,  stagnation  point),  and  simplified  ablation  (i.e., 
pseudo-binary  mixture  of  air  and  ablated  gas). 

One  of  the  earlier  studies  which  accounted  for  the  influence  of  electronic  heat 
conduction  at  high  temperatures  was  that  of  Adams ^  in  which  the  frozen  air 
boundary  layer  with  a  fully  catalytic  wall  was  treated.  Adams’ results, based  on  sim¬ 
plified  transport  properties,  indicated  that  at  entry  apeeds  of  45,  000  ft/sec, 
tlie  convective  heat  transfer  rate  with  ionization  effi>cts  ivas  30  percent  higher 
than  convective  heating  without  ionization  obtained  by  extrapolation  of  existing 
theories. 

To  illustrate  the  significance  of  entry  speeds,  lunar  probes  reentry  the  Earth's 
atmosphere  at  35.  000  ft/sec  while  Mars  probes  reenter  at  45,  000  to  65,  000  ft/sec. 
Air  begins  to  ionize  significamtly  at  less  than  10,  000*K  and  is  completely  single- 
ionized  at  20,  OOO'K.  At  50,  000  ft/sec  and  altitude  of  190,  000  feet,  the  stagna¬ 
tion  temperature  is  15,  000*K. 

Van  Tassell  and  PalloneZ^  employing  similarity-type  solutions,  used  Hansen' 
propi  rtisB  to  study  the  heat  transfer  rate  for  air  in  equilibrium  dissociation  and 
ionization.  Comparing  their  results  and  those  of  Adams  with  Avco  RAD  experi¬ 
mental  results  they  found  good  agreement  between  theory  and  experiment  for 
flight  speeds  up  to  35,  000  ft/sec.  Cohen ■^,  also  using  similarity-type  solutions 
and  employing  Hansen's  properties,  was  able  to  show  cigreement  with  Van  Tassell 
lud  Pallone  for  flight  speeds  up  to  40,  000  ft/sec. 

A  severe  departure  from  the  above-mentioned  theory  was  obtained  by  Scale 
His  results  for  equilibrium  showed  extremely  high  beat  transfer  rates.  The 
iransport  properties  used  by  Scala  were  based  on  a  charge-induced  dipole  model 
for  the  species  N^,  N,  N+,  and  e“.  Pal’one  and  Van  Tassell^  compared  trans¬ 
port  properties  computed  from  Yos^  collision  integrals  based  on  the  species, 

N,  N+,  and  e'  and  were  able  to  show  discrepancies  in  Scala'p  high 

temperature  model.  These  discrepancies  are  discussed  in  detail  in  reference 
6  where  theoretical  predictions  of  thermal  conductivity  are  compared  with  esti¬ 
mates  of  thermal  conductivity  derived  from  liaaecker  ar'  experiments. 
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Ho«hizaki  ,  using  similarity-type  solutions  accounted  for  the  effects  of  disso¬ 
ciation  and  ionization  through  the  useful  concent  of  total  thermodynamic  and 
transport  property  |i.  e. ,  reaction  conductivity,  first  discussed  by  Hirshfelder, 
and  subsequently  used  in  reference  2).  Hoshizaki  employed  Hansen's  properties 
and  obtained  correlations  applicable  to  an  Earth  or  Venus  atmosphere  which 
agreed  well  with  his  own  experimental  heat  transfer  data.  Comparison  of  the 
Hoshizaki  correlation  equation  (equation  20  of  reference  8)  shows  ::^oad  agree¬ 
ment  with  the  results  of  reference  2,  4, and  6,  for  flight  speeds  up  to  40,  000 
ft/ sac.  Adams'  results  for  a  frozen  boundary  layer  were  also  in  good  agree¬ 
ment  with  HoshizakL 

Fay  and  Kemp**  proposed  a  simplified  binary  diffusion  model  in  which  one 
diffusion  equation  as  well  as  three  overall  conservation  equations  are  used  just 
as  is  the  case  of  the  dissociated  gas  boundary  layer.  This  model  is  based  upon 
an  assumption  that  the  relative  velocity  between  the  ions  and  atoms  is  negligible 
compared  with  other  diffusion  velocities.  The  justification  proposed  for  this 
assumption  is  that  the  charge  exchange  cross  section  for  N  -I'T*'  collisions 
(based  on  estimates  by  Yos')  is  much  greater  than  that  for  other  heavy-particle 
collisions,  thereby  impeding  the  relative  motion  of  the  ions  with  respect  to  the 
atoms.  The  results  of  Fay  and  Kemp  for  equilibrium  N2  are  in  good  agreement 
with  those  of  P^allone  and  Van  Tassell  over  a  large  velocity  range. 

The  comprehensive  experimental  results  of  Rose  and  Stankevics^^  using  the 
electrical  driven  shock  tube  confirmed  the  conclusions  of  Pallone  and  Van 
Tassell  and  the  theoretical  results  of  Fay  and  Kemp  for  equilibrixjm  N2  in  addi¬ 
tion  to  verifying  the  theories  of  Cohen  and  Hoshizaki  for  equilibrium  air.  The 
agreement  between  theory  and  experiment  was  thus  fairly  well  established  for 
flight  speeds  to  50,  000  ft/  sec. 

Attaining  equivalent  flight  s-^eeds  substantisdly  higher  than  50,  000  ft/ sec  in  the 
shock  tube  has  proven  difficult  00  that  sufficient  eacperimexxtal  verification  of 
boat  transfer  rates  is  generally  not  available  for  the  very  high  speed  flight 
regime.  Fur&ermore,  uncertainties  in  calciilating  fiie  total  thermal  conductiv¬ 
ity  for  partially  ioidzed  air  have  discouraged  previous  effirts  to  eictend  theoret¬ 
ical  predictions  of  stagnation  point  heat  transfer  to  flight  speeds  above  50,  000 
ft/ sec. 

Howe  and  Schaeffer  have  examined  the  effects  of  uncertainties  in  the  total 
fiierr'tal  conductivity  of  air  on  convective  heat  transfer  for  stagnation  teinpera- 
ttxre  dp  to  30.  000* K.  Using  a  merged  shock -layer  type  of  solution  they  em¬ 
ployed  the  following  variations  of  Yos'  7,  transport  properties:  (1)  Yos' 
thermal  conductivity  of  raference  7,  (2)  Yos'  thermal  conductivity  increased 
by  an  order  of  magnitude  at  the  high  temperature  end.  a»3  (3)  Yos  thermal 
conductivity  faired  into  Hansen's  thermal  conductivity  at  1 3,  000*K.  Their 
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results  indicate  that  the  convective  heat  transfer  rate  is  relatively  insensitive 
to  large  uncertainties  in  the  total  thermal  conductivity.  For  example,  an 
uncertainty  of  a  factor  of  10  in  the  total  thermal  conductivity  of  air  influences 
the  convective  heating  rate  by  only  a  factor  of  1.75  at  a  flight  velocity  of  70,  000 
ft/ sec,  and  by  a  factor  of  2  at  85,  000  ft/ sec.  This  effect  has  also  been  shown 
in  reference  6,  where  a  comparison  wait  made  of  the  heat  transfer  rate  using 
Maecker's  arc  dat^*  and  Yos*  theoretical  data. 

The  object  of  the  present  report  is  to  extend  the  equilibrium  air  analytical 
study  of  Pallone  and  Van  Tassell^  to  flight  speeds  up  to  70,  000  ft/ sec.  For  the 
present  study  we  employ  the  most  recent  transport  and  thermodynamic  prop* 
arties  calculated  by  Yos^^,  and  examine  the  stagnation  point  convective  heat 
transfer.  Results  are  then  compared  with  experiment,  including  recent  Avco 
RAD  shock  tube  data  ^3.  Boundary  layer  profiles  fci  both  two-dimensional  and 
axisymmetric  bodies  are  included  in  this  report.  Results  are  also  presented 
showing  first  order  effects  of  simulated  ablation,  i.  e. ,  pseudo-binary  mixture 
of  air  and  injected  air. 

The  effects  of  blowing  on  the  convective  heat  transfer  has  been  the  subject  of 
many  studies,  including  those  of  Reshotko  and  Cohen^^  who  used  constant- 
property  boundary  layer  solutions.  Libby^^  relaxed  these  restrictions  on  the 
thermodynamic  and  transport  properties  by  taking  the  properties  as  simple 
functions  of  the  temperature.  The  results  of  both  studies,  however,  are  ap¬ 
plicable  only  to  very  'ow  flight  speeds  because  of  the  choice  of  transport  prop¬ 
erties,  Of  special  note  is  an  excellent  report  by  Howe  and  Schaefferl^  which 
e.  iploys  variable  transport  and  thermodynamic  properties  with  a  merged-layer- 
type  solution  to  study  stagnation  point  flows  for  flight  velocities  between  30,  000 
and  50,  000  ft/ sec.  Coupled  radiative  and  convective  heat  transfer  is  presented 
as  wall  as  a  comprehensive  study  of  the  effects  of  moderate  rates  of  injection. 
In  the  present  paper,  the  effects  of  blowing  to  show  first-order  effects  of  abla¬ 
tion  are  examined  at  flight  speeds  up  to  65,  000  ft/  sec  by  extending  the  calcula¬ 
tions  for  the  convective  heat  transfer  rate  to  include  large  rates  of  air  injection 
at  the  stagnation  point.  Both  two-dimensional  and  axisymmetrlc  bodies  are 
studied. 


U.  THERMODYNAMIC  AND  TRANSPORT  PROPERTIES 


Thermodynamic  and  transport  properties  calculated  by  Yos  and  used  in  the 
analysis  of  Pallone  and  Van  Tassell^  have  been  superseded  by  Yos  These 
new  properties,  which  are  employed  in  the  present  report,  are  shown  in 
figures  1  to  8. 

The  new  thermodynamic  properties  are  very  slightly  changed  as  a  result  of 
Yos*  new  calculations.  For  example,  the  specific  heat  for  air  is  only  5  to  15 
percent  higher  at  temperatures  above  15,  000*K.  A  comparison  of  Yos*  new 
thermodynamic  properties  with  those  presented  in  A£DC-TR-64>183  (September 
1964)  showed  excellent  agreement  for  temperatures  up  iu  15,  00C*K,  which  is  as 
high  as  the  AEDC  tables  go. 

Significant  changes  are  indicated,  however,  in  the  transport  properties  of  figures  1 
to  6,  especially  for  the  region  of  partial  ionissation.  At  14,  000*K,  for  example, 
the  total  thermal  conductivity  for  nitrogen  is  approximately  50  percent  higher 
than  Yos*  earlier  predictions  in  reference  7. 

Theoretical  calculations  of  transport  properties  are  based  on  kinetic  theory. 

That  is.  explicit  formulas  for  viscosity  and  thermal  conductivity  as  functions 
of  collision  integral,  particle  mass,  and  mol  fraction  of  particular  species 
are  obtained  from  the  first  Chapman-Enskog  approximation,  or  slight  variation 
thereof  (eed  Yos  ^). 

For  temperatures  up  to  10, 000*K.  the  transport  properties  can  be  predicted 
with  fairly  good  accuracy  so  that  reasonable  agreement  exists  among  various 
investigators.  However,  above  10,000*K.  where  par. ial  ionisation  exists, 
there  has  been  considerable  disagreement  in  the  prediction  of  the  total  thermal 
coudictivlty.  Total  thermal  conductivity,  which  appears  in  the  definition  of  the 
Prandtl  number,  Pr  s  pCpA,  includes  the  effects  of  ordinary  translational 
energy  transport,  thermal  diffusion,  and  chemical  and  ionization  reactions. 

In  the  region  of  partial  ionization  (9.  000  to  16,  000*K)  resonant  charge 'exchange 
processes  in  which  an  electron  is  transferred  from  an  atom  to  a  positive  ion 
of  the  same  species  (e.  g. ,  N  4  N*^-*  N*^  4  N)are  a  major  component  of  the 
thermal  conductivity  of  the  gas.  By  including  the  effects  of  chargo-exchange 
Yos*^  obtained  a  diffusion  cross  section  between  atoms  and  atomic  ions  which 
was  almost  an  order  of  magnitude  larger  than  the  gas'kinetic  cross  section 
used  in  most  previous  high  temperature  transport  property  calculations. 

The  most  recent  calculations  by  Yos^^  yield  charge -exchange  cross  sections 
for  oxygen  which  are  in  excellent  agreement  with  the  experimental  measure¬ 
ments  of  Stebbings,  et  al.  12  experimental  cross  sections  have  yet  been 
measured  for  nitrogen  but  Yos*  recent  theoretical  cross  sections  for  nitrogen 
appear  to  be  in  fair  agreement  with  the  recent  semiemperioal  estimates  of 
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Knof,  Mason,  and  Vanderslice  12,  Yos*  new  nitrogen  cross  section  resiuts  in 
a  calculated  thermal  conductivity  which  is  about  50  percent  higher  than  his 
previous  calculations,  at  the  regions  uf  interest  here.  This  new  prediction  of 
thermal  conductivity  for  nitrogen  is  also  in  good  agreement  with  that  of  Ahtye 
for  temperatures  up  to  15,  000*K.  (See  reference  18.) 

Experimental  verification  of  thermal  conductivity  can  be  made  using  two  gen¬ 
eral  techniques.  The  first,  and  most  commonly  used  technique,  employs  the 
constricted  laminar  arc-column  originated  by  Maeckeri9.  Using  the  constricted 
arc  Maecker  was  able  to  determine  total  thermal  conductivity-  for  nitrogen  at 
atmospheric  pressure  and  temperatures  from  5000  to  15,  OOO’K.  The  second 
experimental  technique  employs  overall  heat  transfer  measurement  in  a  shock 
tube  from  which  one  then  infers  the  values  of  the  total  thermal  conductivity. 

With  the  first  technique,  the  Elenbaas -Heller  equation  for  radial  energy  trans¬ 
port  in  an  optically -thin  cylindrical  arc  column 


is  integrated  once  and  solved  for  the  thermal  conductivity  K,  at  a  rad-  d  dis¬ 
tance  t'  , 


where  £  is  the  electric  field  strength  in  the  column,  r  is  the  radius  of  the  con¬ 
stricting  tube,  9  is  the  electrical  conductivity  of  the  gas,  and  is  the 
power  radiated  from  the  gas  per  unit  volume.  The  gas  properties  a,  K ,  and 

are  generally  functions  of  the  gas  temperature  T,  as  v,ell  as  of  the  nature 
of  tne  gas  and  the  pressure.  The  general  procedure  of  the  measurement  is  to 
first  determine  the  electrical  conductivity  9  (T)  from  the  integrated  form  of 
Ohm's  law,  using  measured  values  of  E,  I  ,  and  the  radial  temperature  dis¬ 
tribution  tC?)  for  several  different  arc  currents.  The  temperature  profiles 
are  obtained  from  continuum  and  line  radiation.  In  the  final  step,  the  thermal 
conductivity  K(T)  !<i  obtained  from  equation  (2)  using  the  previously  determined 
value  of  9  (T)  and  the  measured  value  of  P^^  .  Further  details  of  this  experi¬ 
mental  technique  are  given  in  reference  12. 

Additional  experiments  to  determine  the  thermal  conductivity  of  gases  at 
elevated  temperatures  using  the  constricted  arc  column  have  been  carried  out 
at  Avco  RAD  by  Bennett,  Yes,  Knopp,  Morris,  and  Bade^. 
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Initial  arc  data  obtained  by  the  above  researchers  confirmed  J^4aecker's 
experimental  values  for  nitrogen  at  5000 to  15,  OOO’K.  This  experimertal  arc 
data  when  used  with  the  EUenbaas -Heller  energy  balance  (equation  (1),  ’-alid 
for  optically  thin  gas)  resulted  in  values  of  thermal  conductivity  at  the  higher 
temperatures  which  were  greater  by  an  order  of  magnitude  than  Yos^^  pre¬ 
dictions.  However,  the  experimental  values  were  open  to  question  since  the 
experimental  thermal  conductivity  is  an  apparent  function  of  arc  current. 

This  observation  suggested  to  Bennett,  et  al,  that  some  energy  transport 
mechanism  might  be  present  which  had  not  been  allowed  for  in  the  reduction  of 
the  data. 

A  systematic  analysis  of  the  data  and  careful  research  on  the  abso'  ion  coeffi¬ 
cient  of  argon  and  nitrogen  in  the  vacuum  ultraviolet  led  to  the  conclusion  that 
the  arc  column  emits  significant  radiation  in  the  far  ultraviolet  and  is  relatively 
opaque  to  such  radiation  for  pressures  on  the  order  of  one  atmosphere.  This 
finding  suggested  that  the  assumption  of  an  optically -thin  gas  column  is  not 
applicable  for  use  in  the  £-H  energy  balance. 

When  the  vacuum  ultraviolet  radiation  was  included  in  an  approximate  manner 
in  the  assumed  radiative  loss  for  the  analysis  with  nitrogen,  the  dependence  of 
thermal  conductivity  on  arc  current  was  eliminated.  Furthermore,  the  experi¬ 
mental  values  were  now  found  to  be  in  good  agreement  wi^h  Yos'  predictions, 
as  can  be  seen  from  figure  1.  However,  before  one  can  araw  final  conclusions 
the  eiq>erimeutai  thermal  conductivity  should  be  obtained  from  solutions  of  the 
£-H  equation  with  nongray  radiative  transfer  terms. 

To  get  a  quantitative  understanding  of  the  effect  of  thermal  conductivity  on  the 
heat  transfer  rate,  boundary  layer  solutions  were  obtained  using  the  four 
thermal  conductivity  curves  shown  in  figure  9,  For  each  boundary  layer 
solution  the  values  for  /i  and  were  taken  from  figures  2  and  3,  respectively. 

Figure  10  shows  dramatically  the  effect  of  thermal  conductivity  on  the  heat 
transfer  rate.  In  comparL^g  the  temperature  profiles  for  the  above  cases  it 
is  interesting  to  note  that  the  heat  transfer  rate  is  controlled  primarily  by  con¬ 
ductivity  in  the  vicinity  of  the  wall.  The  dip  in  the  thermal  conductivity  for  the 
bimperature  region  7000  to9000*K  produces  the  effect  of  a  relatively  insulating 
layer  which  divides  the  boundary  layer  into  two  regions,  a  high  temperature 
outer  layer,  and  a  lower  temperature  inner JU'yer.  Values  of  beat  transfer  rate 
tabulated  in  figure  10  also  give  evidence.-that  this  apparent  layer  of  relatively 
low  conductivity  gas  effectively  insulates  the  body  surface  from  the  high  tem¬ 
perature  gases  which  exist  at  the  edge  of  the  boundary  layer. 
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m.  ANALYSIS 


The  well-known  boundary  layer  equations  describing  the  two-dimensional  and 
axisymmetric  flow  of  a  gas  in  equilibrium  dissociation  and  ionization  are  given 
below.  The  effect  of  diffusion  is  contained  in  the  reaction  conductivity  which, 
in  turn,  is  contained  in  the  Prandtl  number.  Transverse  curvature  and  pres¬ 
sure  diffusion  effects  have  been  neglected. 


Continuity  equation 

d  ;  d  . 

— —  (P*r0  u  0 

o*  dy 


(3) 


Momentxnn  equation  in  x-direction 


da  da 

P“  -r-  +  rr- 

dx  ay 


£p 

dr 


Energy  equation 


dH 

oa  +  pa 
dx 


dH 


d(u^/2)  I 

dy  I 


(4) 


(5) 


Input  tables  of  thermodynamic,  properties  are  employed  in  the  present  analysis 
in  place  of  the  equation  of  state.  In  addition,  the  analysis  of  the  injection 
region  can  be  made  of  practical  interest  by  adding  the  heat  balance  relation, 


Introducing  the  Levy-Lees  trant>formation  and  the  concept  of  local  similarity, 
the  above  equations  reduce  to  the  following  form; 


(7) 

(8) 
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where. 


The  appropriate  bouadary 
At  1}  •>  0  :  f  M  f  (0)  “ 

At  tf  ■*  m:  -•  1  ; 


Pm  M*  «i 


conditions  for  these  equations  are: 

- ^  ;  £,-0 

PmPm  »I 


i  ■*  i 


In  addition,  if  the  energy  balance  is  required  as  an  optioital  boundary  condition, 
the  following  equation  holds. 


*w“«c 


%f(0) 


(9) 


If  there  is  no  injection,  or  if  tite  use  of  an  energy  balance  is  not  desired,  then 
either  or  [gj  is  given  and  the  remaining  boundary  values  are  found  in  the 
solution^  * 


DI^CRIPTION  OF  PKOGlkAM 


Particular  solutions  to  equations  (7),  (8)  and  (9)  were  obtained  by  employing 
the  iterative  integration  metttod  described  in  detail  by  Van  TasscU  and  Pallone 
In  this  scheme,  an  additional  transformation  is  introduced  which  reduces  the 
above  system  of  ordinary  differential  equations  to  a  set  of  five  first-order  dif- 
ferentiai  equatione  which  are  then  solved  by  an  Iteraticm.  procedure  employing 
1000  points  across  the  boundary  layer.  The  Prandtl  number  and  pg  ratio  are 
allowed  to  vary  continuously,  and  convergence  is  very  rapid.  A  typical  boimdary 
layer  solution  requires  approximately  one-half  minute  of  machine  time  with 
usually  less  than  10  iterations  per  case. 
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IV.  RESULTS  AND  DISCUSSION 

Boundajry  layer  characteristics  (including  velocity  and  temperature  profiles, 
and  heat  transfer  rates)  were  computed  at  the  stagnation  point  for  a  spherical 
and  cylindrical  body  with  a  radius  of  one  foot  and  a  representative  stagnation 
pressure  of  one  atmosphere  and  10  atmospheres.  The  calculations  covered  the 
range  of  flight  speeds  from  5000  to  70,000  ft/sec.  Values  of  the  blowing 
parameter.  -  ,  included  0.  0.  0.25,  0.50,  and  1.0.  The  injected  fluid  was 

assumed  to  be  equilibrium  air. 

Table  1  summarizes  the  results  for  the  cases  studied  in  this  report.  The  bound¬ 
ary  layer  profiles  for  these  cases  are  given  in  tables  II  through  VII.  For 
laminar  boundary  layers  it  is  well  known  that  q  varies  with  die  square  root  of  the 
pressure  .  Since  die  transport  properties  are  known  to  exhibit  a  very  weak 
dependence  on  pressure,  it  follows  that  the  heat  transfer  rates  of  table  1  can  be 
scaled  for  arbitary  stagnation  pressures.  Typically,  such  scaling  would  be  valid 
for  a  Mars  reentry  trajectory  where  stagnation  pressures  may  vary  from  0.01 
atmosphere  to  several  atmospheres. 

It  may  be  noted  in  table  I  that  when  blowing  is  present,  and  in  particular  if 
is  kept  constant,  then  q  as  a  function  of  flight  speed  wilt  go  through  a  maximum. 
This  peculiar  behavior  of  4  should  not  be  interpreted  in  the  practical  sense 
that  4  will  decrease  for  constant  blowing  with  increasing  V„.  Rather,  this 
behavior  in  4  i*  &  result  of  our  definition  of  and  of  the  constraint  that 
remain  constant  with  increasing  If  one  defines  the  blowing  parameter  as 
fj  =  -  so  that  Pit  is  evaluated  at  conditions  at  the  edge 

of  the  boundary  layer,  rather  than  at  the  wall  as  is  done  in  this  report,  then 
4  will  increase  monotonically  with  increase  in 

The  calculated  heat  transfer  rates  for  zero  blowing  are  presented  in  figure  1 1 
for  comparison  with  recent  shock  tube  measurements.  Since  we  are  interested 
mainly  in  the  flight  regime  above  40,  000  it/ see  the  numerous  theoretical  results 
cited  in  the  literature  which  are  applicable  to  flight  speeds  below  40,  000  ft/sec 
have  not  been  included.  It  may  be  noted,  however,  that  Hoshizaki’s  correla¬ 
tion  equation  for  cquilibritnn  air  (equation  20  of  reference  8)  based  on  Hansen's 
transport  properties  would  lie  somewhat  below  our  curve  at  =  20,  000  ft/ sec 
and  would  intersect  our  curve  at  =  44,  000  ft/sec. 

For  clarity,  only  the  curve  fits  from  Rose  auad  Stankevicsl®  are  shown  in 
figure  11  although  their  curve  fits  represent  a  considerable  number  of  data 
points  extending  to  equivalent  flight  speeds  of  50,  000  ft/ sec. 

The  Avco  RAD  shock  tube  data  shown  in  figure  11  indicates  that  shock  velocities 
very  near  13  mm/  gsec  (equivalent  flight  velocity  of  56,  000  ft/sec)  are  now 
being  attained  in  the  electrical  driven  shock  tube.  At  the  flight  speeds  above 
50, 000  ft/sec.  the  Avco  RAD  experimental  data  falls  slightly  above  our 
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calculated  results.  Because  the  test  models  have  a  relatively  small  nose 
radius  it  seems  unlikely-  that  radiation  from  the  surrounding  gas  cap  could  be 
significant  at  the  prevailing  test  conditions.  It  is  possible,  however,  that  the 
model  may  be  seeing  ultraviolet  radiation  from  the  slug  of  test  gas.  Experi¬ 
ments  are  currently  in  progress  to  accurately  determine  the  relative  contri¬ 
butions  due  to  radiation  from  the  slug  of  test  gas  and  from  the  surrounding 
gas  cap. 

The  stagnation  point  velocity  gradient  used  in  the  calculation  of  heat  transfer 
rate  for  the  two-dimensional  bodies  was  based  on  Newtonian  theory^^  and  was 
calculated  to  be  approximately  6  percent  larger  than  the  stagnation  point  velocity 
gradient  for  an  axisymmetric  body.  This  value  for  the  two-dimensional  body 
has  been  verified  experimentally  by  Korkan^l^ 

Theoretical  predictions  of  the  convective  heat  transfer  rate  for  equilibrium 
nitrogen  are  shown  in  figure  12  and  compared  with  recent  Avco  RAD  shock 
tube  data.  The  agreement  is  excellent  at  flight  speeds  between  30,  000  and 
40, 000  ft/sec.  Not  enough  data  is  available  at  the  flight  speeds  above  50,  000 
ft/sec  to  clearly  show  the  trends  indicated  in  the  experimental  data  of  figure  11. 
Theoretical  heat  transfer  rates  calculated  by  Fay  and  Kemp^  using  a  simplified 
binary  diffusion  model  are  shown  in  figure  12  to  be  only  slightly  higher  than 
those  obtained  in  the  presen*  eport.  It  can  still  be  said,  however,  that  in 
figure  12  the  agreement  bet’  .n  both  theories  and  experiment  is  very  good. 

Results  obtained  for  axisymmetric  bodies  using  a  pseudo-binary  mixture  of 
air  and  injected  air  are  plotted  in  figure  13.  The  effect  of  the  blowing  parameter 
is  seen  to  be  seen  to  be  strongly  nonlinear  at  the  higher  blowing  rates. 

In  order  to  obtain  a  quantitative  understanding  of  the  effect  of  transport  prop¬ 
erties  on  the  heat  transfer  rates  comparisons  have  been  made  with  previt^us 
investigators  as  shown  in  figure  14.  The  results  of  Reshotko  and  Cohen^^  for 
a  constant-property  boundary  layer  are  shown  for  two  different  temperature 
ratios.  When  the  constant-property  restriction  is  relaxed  slightly,  as  for 
example,  in  the  result  shown  for  Libby^^  where  the  transport  and  thermo¬ 
dynamic  properties  are  taken  as  simple  factions  of  the  temperature,  then 

the  heat  transfer  parameter  Nu/v^  is  reduced.  When  the  exact  thermodynamic 

and  transport  properties  are  used,  figure  14  shows  that  the  level  of  Su/y/Wc  is 
even  further  reduced,  amd  in  addition,  a  velocity  dependence  is  introduced. 
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Fl««*  12  COMPARISON  OF  THEORETICAU  PREDICTIOHS  WITH  SHOCK  TUBE  DATA  FOR  MITROGEM  WITH  ZERO 

BLOWING 
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\  Rgur*  13  EFFECT  BLOWING  ON  THE  HEAT  TRANSFK  PARAMETER  AT  THE  STAGNATION  POINT 

!  OF  AN  AXISY\WETRIC  BLUNT  BODY  (AIR-AIR) 
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HEAr  TRANSFER  fARAMGTER 
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BLOWmG  MRAMETER,-fw 


Rgor*  :4  COMPARISON  OF  RESULTS  SHOWING  VARIATION  OF  NUSSELT  NUMSBi  WITH  BLOWiNG 
AT  THE  STAGNATION  POINT  OF  AN  AXISYMMETRIC  BLUNT  BODY  (AIR-AIR) 
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V.  CONCLUSIONS 


1.  Discrepancies  between  experiment  and  theory  have  been  ..ignificantly 
reduced  both  by  an  improvement  in  Vos'  predictions  and  aiso  by  a  further 
refinetaent  in  the  data  reduction  from  arc  experiments. 

2.  The  heat  tsansfer  rate  to  the  wall  is  strongly  affected  by  the  shape  and 

the  occurance  of  the  dip  in  the  graph  of  thermal  conductivity  versus  temperature, 

3.  Good  agreement  between  theory  and  experiment  has  been  obtained  for 
equilibrium  air  for  flight  velocities  up  to  S5,  000  ft/sec;  for  equilibrium 
nitrogen  agreement  between  theory  and  experiment  is  excellent  up  to  flight 
velocities  of  50,  000  ft/sec. 
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